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ABSTRACT 

Formulations  for  predicting  dynamic  and  acoustic  responses  of  a  finite  plate  flush-mounted 
on  an  infinite  baffle  subject  to  turbulent  and  mean  flow  excitations  are  presented.  In  deriv¬ 
ing  these  formulations,  the  effects  of  stretching  due  to  in-plane  force  and  coupling  between 
structural  vibration  and  acoustic  radiation  are  taken  into  account.  The  resulting  equation 
governing  plate  vibration  contains  cubic  nonlinearities.  To  obtain  am  approximate  solution, 
the  plate  flexural  displacement  is  expanded  into  orthogonal  base  functions.  The  unknown 
coefficients  associated  with  base  functions  are  determined  by  solving  a  set  of  coupled  nonlin¬ 
ear  integral  equations  using  Galerkin’s  method.  A  stability  analysis  is  given  using  the  basic 
existence-uniqueness  theorem.  In  particular,  stable  conditions  for  a  linearized  system  are 
obtained  using  the  Routh  algorithm.  It  is  shown  that  structural  instabilities  can  be  induced 
by  fluid  loading  and  mean  flow.  Two  unstable  mechanisms  are  found  to  be  attributable  to 
added  damping  and  stiffness  due  to  acoustic  radiation.  The  effect  of  added  damping  in¬ 
creases  linearly,  while  that  of  added  stiffness  increases  quadratically  with  the  Mach  number 
of  mean  flow.  Finally,  the  cross-power  spectral  density  functions  of  plate  flexural  vibration 
and  radiated  acoustic  pressure  are  derived  and  expressed  in  terms  of  the  cross-power  spectral 
density  function  of  turbulent  flow. 


‘This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Con¬ 
tract  No.  NAS1-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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INTRODUCTION 

Responses  of  elastic  structures  subject  to  turbulent  and  mean  flow  excitations  have 
long  been  of  great  practical  engineering  importance,  for  example,  in  the  design  of  aircraft 
and  in  the  reduction  of  aerodynamic  noise  transmission  into  an  aircraft  cabin.  Experi¬ 
mental  and  theoretical  investigations  on  turbulent  flow  and  its  interaction  with  an  elas- 
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tic  structure  have  been  carried  out  by  many  people,  for  example,  Ludwig  ,  Corcos  ’  , 
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Maestrello  ’  ,  Maestrello  and  Linden  ,  Farabee  and  Geib  ,  Martin  and  Leehey  ,  Yen  et 

9  10  11  12 

al.  ,  Ffowcs  Williams  ,  Gedney  and  Leehey  ,  Farabee  and  Casarella  ,  etc.  Most  of 

the  early  work,  however,  has  focused  on  the  behavior  of  linear  systems.  Recently,  with 

the  development  of  modem  computational  facilities,  nonlinearities  in  an  elastic  structure 

and  in  acoustic  wave  propagation  can  be  accounted  for,  and  the  dynamic  response  of 

an  elastic  structure  subject  to  a  turbulent  boundary  layer  excitation  can  be  calculated 

directly^’  Nevertheless,  the  physics  involved  in  interaction  between  turbulent  flows 

and  vibrating  structures  is  still  not  well  understood,  and  usually  cannot  be  revealed  by 

direct  numerical  simulations. 

In  this  paper,  we  derive  an  analytical  formulation  for  calculating  the  dynamic  and 
acoustic  responses  of  a  finite  plate  clamped  flush  to  an  infinite  baffle  under  the  excitation 
of  turbulent  flow,  with  the  structural  nonlinearities  due  to  stretching  of  in-plane  force 
taken  into  account  (see  Section  I).  Because  of  the  presence  of  nonlinearities  and  coupling 
of  structural  modes  to  the  radiated  acoustic  pressure,  closed  form  solutions  cannot  be 
found.  In  order  to  obtain  an  approximate  solution,  the  plate  flexural  displacement  is 
expanded  into  the  orthogonal  base  functions  in  the  plate’s  longitudinal  and  transverse 
directions.  The  unknown  coefficients  associated  with  these  base  functions  are  determined 
by  Galerkin’s  method.  These  are  done  in  Section  II.  In  Section  III,  we  present  an  analysis 
of  structural  instabilities  induced  by  fluid  loading  and  flow.  The  objective  of  this  stability 
analysis  is  to  gain  physical  insight  into  this  complex  problem  by  relating  the  instability 
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phenomenon  to  the  physical  quantities  possessed  by  an  elastic  structure.  For  steady-state 
responses,  we  demonstrate  that  the  cross-power  spectral  density  function  of  the  plate 
flexural  displacement  can  be  expressed  as  the  cross-power  spectral  density  function  of 
turbulent  flow  (see  Section  IV).  Based  on  the  plate  flexural  vibration  response,  we  can 
then  derive  a  formulation  for  calculating  the  cross-power  spectral  density  function  of  the 
radiated  acoustic  pressure  (see  Section  V).  Some  concluding  remarks  are  made  in  Section 
VI. 


I.  PLATE  EQUATION 

Assume  that  a  plate  of  length  L  and  width  b  is  clamped  on  its  edges  to  a  rigid  and 
infinitely  extended  baffle.  On  one  side  of  the  plate,  the  flow  is  turbulent,  with  boundary 
layer  thickness  6.  Outside  this  boundary  layer,  the  fluid  moves  at  mean  velocity  U  (see 
Fig.  1).  The  other  side  of  the  plate  is  assumed  to  be  vacuum,  so  that  the  effect  of  fluid 
loading  acts  on  one  side  of  the  plate  only. 

To  carry  out  the  analysis,  we  make  the  following  assumptions:  (1)  the  amplitude  of 
the  plate  flexural  vibration  is  small,  so  that  the  boundary  layer  will  not  be  significantly 
altered;  (2)  the  turbulent  boundary  layer  is  stationary  in  time  and  homogeneous  in  space, 
so  that  the  pressure  fluctuations  can  be  expressed  as  a  cross-correlation  function  which 
decays  with  spatial  and  temporal  separations  and  convects  downstream  with  the  flow;  and 
(3)  the  pressure  fluctuations  exerted  on  the  plate  surface  due  to  a  turbulent  boundary  layer 
and  due  to  radiated  acoustic  pressure  are  additive. 

Under  these  conditions,  we  can  write  the  equation  governing  the  plate  flexural  vibra¬ 
tion  as 


nf  &  ,  &  V  /  .  L^2w(x,y,t) 

+  dxi  -  +7-  +P,h  w  - 


a®3 

=  -p(x,y,0,t), 


(1) 


2 


subject  to  the  boundary  conditions 


u>(0,  y,  t)  =  w(L,  y,  t)  =  w(x,  0,  t)  =  w(x,  b,  t)  =  0, 

=  0. 

,=6 


dw 

dw 

dw 

dw 

dx 

r=0 

dx 

x=l  dy 

,=o 

(2a) 

(») 


In  Eq.  (1),  the  quantity  D  =  Eh3/ 12(1  -  p2)  is  the  plate  flexural  rigidity,  E  the  Young’s 
modulus,  n  the  Poison’s  ratio  of  the  material,  h  the  thickness,  pp  the  mass  density,  7 
the  structural  viscous  damping  coefficient,  w  the  plate  flexural  displacement  in  the  2-axis 
direction,  and  Nz  the ’stress  resultant  in  the  x-axis  direction  due  to  stretching  of  plate 


bending  motion 
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L  dw(x,y,t ) 
dx 


dx. 


(3) 


The  effects  of  the  stress  resultant  in  the  y-axis  direction  Ny  and  that  in  the  tangential 
direction  Nxy  are  assumed  small  compared  with  Nx ,  and  are  therefore  neglected  in  Eq. 
(1). 


The  term  on  the  right  side  of  Eq.  (1)  depicts  the  excitation  forcing  field  due  to 
turbulent  flow.  Inside  the  turbulent  boundary  layer,  the  quantity  p  can  be  expressed  as 
an  integral  representation/  with  a  volume  integral  representing  the  contribution  from 
distributed  quadrupoles,  plus  a  surface  integral  representing  the  effect  of  sound  reflection 
and  diffraction  on  the  surface.  However,  the  stress  tensor  in  the  volume  integration 
contains  quantities  that  are  unknown  until  the  entire  flow  field  is  solved,  which  is  not 
possible  in  most  engineering  applications.  Hence,  a  different  approach  must  be  taken. 


Note  that  the  effect  of  turbulent  flow  is  to  excite  a  structure  and  generate  sound. 
Therefore  we  can  assume  that  the  pressure  fluctuations  acting  on  the  plate  surface  p  can 
be  expressed  as  a  sum  of  two  parts 


p(z,  y,  O’,  t)  =  pr(x,  y,  0,  t )  +  p(x,  y,  0,  <), 


(4) 
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where  pr  represents  a  forcing  field  which  exerts  the  same  amount  of  pressure  fluctuations 
on  the  structure  as  the  turbulent  boundary  layer,  and  p  depicts  the  effect  of  acoustic 
radiation.  Here  we  assume  that  the  forcing  field  pr  is  random  in  nature,  but  stationary 
in  time  and  homogeneous  in  space.  Therefore  pr  can  be  described  by  a  space-time  cross- 
correlation  function  and  determined  experimentally.  The  function  p  satisfies  the  convected 
wave  equation  outside  the  boundary  layer.  Solution  for  p  can  be  facilitated  by  introducing 
a  potential  function  which  satisfies^ 


1(1 
c2  \dt 


x,y,z,t)  =  0, 


(5) 


subject  to  the  condition  of  continuity  of  normal  components  of  the  fluid  particle  velocities 
at,  strictly  speaking,  the  interface  between  the  turbulent  boundary  layer  and  mean  flow. 
However,  it  is  impossible  to  determine  this  interface  exactly  because  the  flow  beneath  it  is 
turbulent.  To  simplify  the  problem,  we  set  this  interface  at  the  plate’s  nominal  position 
so  that  the  boundary  condition  becomes 


d<j>{x,y,  z,t ) 
dz 


z=0 


(&  +  u£)w(x,y,t)  xeA] 

0  otherwise; 


(6) 


where  A  stands  for  the  plate  surface  area,  U  is  the  mean  flow  speed,  and  c  is  the  speed  of 
sound  of  the  fluid  medium. 


Solution  for  <j>  can  be  obtained  by  taking  a  temporal  Fourier  transformation  of  Eqs. 
(5)  and  (6)  into  the  frequency  domain,  followed  by  a  spatial  Fourier  transformation  into 
the  frequency- wavenumber  domain.  Once  the  frequency- wavemunber  domain  solution  that 

A 

represents  an  outgoing  wave  is  obtained,  the  frequency  domain  solution  of  <j>  can  be  ob¬ 
tained  by  taking  an  inverse  spatial  Fourier  transformation 

fax,V,z,u)  =  “(2^/  J0  G'(x’V»z!z,>y'»°)  (-lu,  +  t7^)  ti(a;',y',w)dx'dy',  (7) 
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where  the  symbol  “  over  a  function,  say,  <j>  implies  a  Fourier  transformation  of  that  function 
defined  by 


^(x,y,t)=  f  4>(x,y,w)eiut<Ju, 

J  — oo 

J  <f>(x,y,t)e~iutdt. 


(8a) 

(86) 


The  quantity  G  in  Eq.  (7)  is  a  Green’s  function.  For  a  supersonic  flow,  the  acoustic 
pressure  is  nonzero  only  when  the  field  point  x  is  inside  the  Mach  cone.  Hence  we  obtain0 


..  *2* cos{tfc[A/(x  -  x')  +  Rauper\/(M 2  -  1)} 

G{x,y,z\x  , y  ,0)  = - - - - 

where  Rsnper  is  the  distance  between  the  source  and  receiver  inside  the  Mach  cone 

R»up  er  =  \/(*  -  -  (Je  -  i)[(y  -  V)1  +  7} 


(9) 


(10) 


Accordingly,  G  for  a  subsonic  flow  can  be  written  as 


G(x,y,z\x',y',0)  = 


fl.«4V,l  -  Af1 


where 


Rtuk  =  V(x  ~  x')2  +  (i  -  A^)[(v  -  y1)2  +  **J 


where  M  is  the  Mach  number  of  the  mean  flow. 


(11) 


(12) 


The  frequency  domain  solution  for  the  radiated  acoustic  pressure  is  then  given  by 


p(z,y,z,u)  =  -/>(-iu  +  V^\  i(x,y,z,u). 


(13) 


The  corresponding  solution  for  the  acoustic  pressure  on  the  plate  surface  in  the  time  domain 
p(x,y,Q,  t)  can  be  obtained  by  taking  an  inverse  Fourier  transformation  of  p. 
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Interchanging  the  order  of  integrations  with  respect  to  the  frequency  and  spatial 
domains,  and  carrying  out  the  integration  with  respect  to  u;  first,  we  obtain 

*  ,L'  (i  +  ve,.,.r)  aMx-y.t1) 


VM*-1  dt’2 


dx'dy' 


t'=r 


+ 


f  fL\ _ h 

Jo  Jo  .Raupery 


T  auper  dw(x',y',t')  (1  +  UQ,uper)  d*w(x',y',t') 


pU  _ 

(2*)2  Jo  Jo  U,«perVM^TI  dt> 


RsuperVM*  -  1  dx'dt1 


dx'dy ' 


t'-r 


PoU2 

(27t)2 


6  fL 


fl 


T,«fer  dw(x',y',t') 
RsuptrVM*  -  1  dx' 


dx'dy', 


(14) 


«'=r 


for  supersonic  flow.  The  integrands  in  Eq.  (14)  are  to  be  evaluated  at  the  retarded  time 


t  =  t  —  Attuptr,  with  At  defined  as 


A t,uper  =  k[M(x  -  x')  +  R,uper]/(M 2  -  1), 


The  quantities  0,.per  and  T..per  in  Eq.  (14)  are  given  by 


(15) 


e 


auper 


c(M 3  - 1) 


(x  -  x') 


R. 


auper 


(x  -  x') 

•l  tuptr  —  ~^2  • 

•***«Per 


Similarly,  for  subsonic  flow  we  have 
r*  rL 


(16a) 

(166) 


f K*’v’°’t)  -  W  L  l  iwi  -A/*  — 5?r — J  „  J1  dy 


£L  r  h_ 

(2*)2  J0  J0  [f2, 

PoU 2  f*  [L 

(2*)2  Jo  Jo 


*  /ir  Ttwt  dw(x',y',i')  (1  +  Ue,ui)  &w(x',y',t') 

by/l  -  M2  dt< 


T,.* 

dw(x',y',t') 

-  M2 

dx ' 

R,nbVl  -  M2  dx'dt' 
dx'dy', 


dx'dy' 


i'=r 


(17) 


t'=r 


where  r  =  t  —  At,«*  with  A t««*  defined  as 

At..*  =  k[M(x  -  x')  +  fJ..*]/(l  -  Af2). 


(18) 
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The  quantities  atub  and  bauf,  in  Eq.  (17)  are  given  by 


1 

”  c(l  —  M2)  . 


(x-  x') 


( x-x ') 


(19a) 

(196) 


Substituting  Eq.  (4)  into  (1),  we  can  rewrite  the  plate  equation  as 

+  £)'•<•.  ».*)  -  + + Pph^l 


=  - pr{x ,  y,  0,  t)  -  p(x,  y,  0,  i),  (20) 

where  pr  and  p  are  fully*  defined. 

II.  GALERKIN’S  METHOD 

Equation  (20)  represents  a  nonlinear,  integral-differential  equation.  Analytical  so¬ 
lutions  to  such  an  equation  cannot  be  found  in  general.  Hence  in  this  section,  we  seek 
approximate  solutions  via  Galerkin’s  method. 

Let  Wk(x)  be  the  kth  eigenfunction  for  a  clamped  plate  in  the  longitudinal  (x-axis) 
direction  and  W<(y)  be  the  zth  eigenfunction  in  the  lateral  (y-axis)  direction  given  by 

Wjfc(x)  =  a*  [sin(Afcx/L)  —  sinh(Afcx/L)]  +  [cos(Atx/L)  —  cosh(Ajtx/L)] ,  (21a) 

W»(y)  =  o’.  [sin(A<y/6)  -  sinh(A<y/6)]  +  [cos(A<y/6)  -  cosh(Ajy/6)] ,  (216) 

where  at  is  the  kth  modal  ratio  defined  as 


<*k  = 


sin  A*  4*  sinh  A* 
cos  A*  —  cosh  A*  ’ 


where  A*  is  the  corresponding  kth  eigenvalue  determined  by 


cos  A*  cosh  A*  =  1. 


For  example,  Aj  =  4.73,  A2  =  7.8532,  and  A3  =  10.9956.  For  a  large  k,  A*  — ►  (2k  +  l)n/2. 


It  is  noted  that  the  eigenfunctions  Wu{x)  are  orthogonal  to  each  other 


£  Wk(x)W,(x)dx  =  |  £  *  -/’  (24) 

However,  the  products  of  PT*(x)  and  its  derivatives  are  not  necessarily  orthogonal. 
Similar  relationships  also  hold  true  for  Wi(y) 

jy,(y)W,(y)iy  =  {l  \=j}  (25) 

Next,  we  expand  w(x,y,t)  in  terms  of  the  orthogonal  base  functions  (W(x)}  and 
{W’(y)}  which  satisfy  the  boundary  conditions,  Eq.  (2), 

»(*,¥,<)  =  W*)}7'[C(t)HW'M}  (26) 

where  (W^x)}  and  {P^(y)}  represent  column  vectors  containing  the  normal  modes  Wk(x) 
and  Wi(y),  respectively.  A  superscript  T  in  Eq.  (26)  indicates  a  transposition,  and  [C(t)] 
represents  the  matrix  consisting  of  the  unknown  coefficients  associated  with  {W(x)}  and 
{iy(j/)}.  Physically,  the  jth  column  vector  {C'7}  represents  the  amplitude  of  coupling 
between  all  the  normal  modes  in  the  x-axis  direction  and  the  jth  mode  in  the  t/-axis 
direction. 

To  facilitate  the  derivation  of  solution  formulation,  we  also  expand  pr(x,  y,  0,  t),  which 
represents  the  forcing  field  due  to  a  turbulent  boundary  layer,  by  (VT(x)}  and  {W"(y)} 

w(»,  Vi  o.  t)  =  {W'(x)}7'(/(()]{IV(»))  (27) 

where  [f(t)]  represents  the  amplitude  of  the  excitation  forcing  field. 

Substitute  Eqs.  (26)  and  (27)  into  (20),  multiply  the  resulting  equation  by  |W(x)}, 
integrate  over  x  from  0  to  X,  and  then  multiply  the  equation  by  {W(y)}T  and  integrate 
over  y  from  0  to  b.  This  procedure  then  leads  to  the  following  matrix  equation 

MO]  +  [*][C]  +  [X][C)  =  [/],  (28) 
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where  [$],  [#],  and  [x]  consist  of  sub-matrices  [$,J],  [¥*■*],  and  respectively,  which 
represent  the  effects  of  coupling  between  all  the  modes  in  the  x-axis  direction  and  the  jth 
mode  in  the  y-axis  direction.  The  elements  of  each  of  the  sub-matrices  [$'•>],  [x,jf], 

and  [S'y]  can  be  written  in  the  following  general  forms 


$*/  =  PphSijSki 


+iV$bL  /  Jo  il  Jo  ^  ^  dx'dy'dxdy  (29.) 

l\vk(x')  -  (1  WMWMW^x)  dx'dy'dxdy  (29b) 

Xg  =  D{[(t)<+  (t)4]  6‘ih‘  +  W^WMdy  £  Wk(x)W,(x)dz\ 

P"U2  £  £  £  \^g~-Wi(y')Wj(y)Wl(x)dx'dy'dxdy  (29c) 


•="J  _ 

“w  2  bL2 


(2n)HL 

jgfo  rb  rt  I  rl<  T  TO  n 


L  L  \l 


2 

dx' 


x Wi(») w> (»)  ", WK*)** 


(29d) 


where  and  are  Kronecker  deltas.  For  a  supersonic  flow,  we  can  simply  replace 
quantities  R,  9,  and  T  in  Eq.  (29)  by  RaUper >  0*uper,  and  T,ttper  given  in  Eqs.  (10)  and 
(16),  respectively.  Similarly,  for  a  subsonic  flow,  we  can  replace  these  quantities  by  Rsub, 
Stub,  and  Ttub  given  by  Eqs.  (12)  and  (19),  respectively. 

Physically,  the  matrices  [$],  [$],  [x],  and  [E]  in  Eq.  (28)  reflect  the  effects  of  mass, 
damping,  and  stiffness  per  unit  area  of  the  plate,  respectively.  In  particular,  the  quadruple 
integrals  involved  in  the  elements  of  the  sub-matrices  and  x*j  represent  the  effects 

of  added  mass,  added  damping,  and  added  stiffness  per  unit  area  due  to  the  radiated 
acoustic  pressure  field,  respectively.  The  elements  E^  are  nonlinear  and  contain  quadratic 
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powers  of  the  unknown  coefficients  {C,J}.  Note  that  in  this  case  coupling  occurs  not  only 
through  nonlinearities,  but  also  through  added  mass,  added  damping,  and  added  stiffness 
due  to  acoustic  radiation. 

The  elements  of  the  sub-matrices  [/*■*]  on  the  right  side  of  Eq.  (21)  are  given  by 

fli = -Skim-  (30) 

Note  that  in  order  to  simplify  the  problem,  we  have  assumed  that  the  time  delay  At 
is  negligibly  small  in  the  derivation  of  Eq.  (28).  This  assumption  is  permissible  if  the  plate 
has  a  finite  size  and  if  the  time  required  for  an  acoustic  signal  to  traverse  the  plate  surface 
is  substantially  less  than  the  typical  period  of  the  signal. 

III.  STABILITY  ANALYSIS 

Equation  (28)  represents  a  system  of  coupled,  nonlinear,  integral  equations.  Because 
of  the  presence  of  nonlinearities  and  coupling  between  structural  vibration  and  acoustic 
radiation,  the  amplitude  of  the  plate  flexural  vibration  may  become  unstable  under  certain 
circumstances.  Therefore  a  general  stability  analysis  is  given  below. 

Equation  (28)  shows  that  for  N  longitudinal  modes  and  M  transverse  modes,  one 
need  solve  an  ( N  x  M)  x  (N  x  Af)  matrix  equation.  Solutions  to  such  a  matrix  equation 
are  extremely  involved.  For  the  purpose  of  illustration,  let  us  consider  only  one  mode  in 
the  y-axis  direction  and  two  modes  in  the  x-axis  direction.  Thus,  Eq.  (28)  reduces  to  a 
2x2  matrix  equation.  The  stability  of  such  a  system  can  be  analyzed  by  examining  the 
homogeneous  part  of  that  equation 

me}  +  [*]{£}  +  lx){C}  -  [E ){C]  =  0  (31) 


10 


where  the  elements  of  matrices  [$],  [¥],  [x],  and  [E]  are  given  by  Eqs.  (29),  with  the 
superscripts  omitted  for  brevity.  In  particular,  the  elements  of  the  matrix  [E]  are  quadratic 
in  C  and  can  be  written  as 

=  (efc/Ci  +  CkiCi)2  (32) 

where  the  elements  of  and  Cki  are  defined  as 

HI  =  £  W‘'(x)W,(x)ix  J°L  W?(x')dx'  w‘(y)<{y  (33a) 

(kl  =  W''(x)W,(x)dx  £  W?(z)dx  jf  ‘  W;(y)iy  (33 h) 

with  indices  k  and  /  varying  from  1  to  2. 

19 

Now  we  use  the  basic  existence-uniqueness  theorem  to  analyze  the  stability  of  the 
system  defined  by  Eq.  (31).  First,  let  us  rewrite  Eq.  (31)  in  the  following  way 


{C}  =  [a){6}  +  [fi){C)  +  [u){C} 


where  matrices  [a],  [/?],  and  [u]  are  defined  as 

M  =  -  WM*] 

m  =  -m-'ix] 

M  =  W_1[H] 


where  [$]  1  represents  the  inversion  of  the  mass  matrix. 
Next,  we  define  new  variables 


Y1=Cl 
Y2  =  Ci 

y,  =  c2 

n  =  c2. 


(34) 


(35a) 

(356) 

(35c) 


(зба) 

(збб) 
(36c) 
(36d) 
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Substituting  Eq.  (36)  into  (34),  we  obtain 


{Y)  =  {J[A]  +  J[B]){Y} 


(37) 


where  (Y)  is  the  time  derivative  of  {Y}  given  by  {Y}T  =  {Yj,  Yj,  Y3,  Y4},  J[A\  is  the 
Jacobian  of  the  linear  part  of  Eq.  (34)  defined  by 


J[A  = 


■  0 

1 

0 

0  - 

Ai 

an 

Al2 

an 

0 

0 

0 

1 

u  An 

0(21 

P22 

022- 

(38) 


where  a*j  and  /?*!  are  given  by  Eq.  (35).  Similarly,  J[B]  in  Eq.  (37)  is  the  Jacobian  of 
the  nonlinear  part  of  Eq.  (34)  defined  by 


J[B)  = 


0  0  0  on 

(«M*i+Cntt)*  0  (elaY,  +  C12Y3)2  0 

0  0  0  0 

L(€2iY,+C2in)2  0  (enYi+CnY*)2  Oj 


(39) 


where  e*j  and  C ki  are  given  by  Eq.  (33). 

The  stabilities  of  the  system  defined  by  Eq.  (37)  can  now  be  examined  by  solving 


for  the  equilibria  obtained  by  setting  the  left  side  of  Eq.  (37)  to  zero 

(J[A)  +  J[B)){Y}=0. 
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(40) 


In  general,  for  a  nonlinear  system  with  cubic  nonlinearities  there  are  multiple  equilibrium 
20 

positions,  and  the  plate  flexural  vibration  may  become  chaotic  self-excited  oscillations 
in  the  presence  of  mean  flow. 

For  simplicity,  we  demonstrate  in  what  follows  a  stability  analysis  of  a  linearized 
system  obtained  by  omitting  the  nonlinearities  in  Eq.  (37) 

{Y}  =  J[A){Y).  (41) 
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Without  loss  of  generality,  we  can  assume  a  form  of  solution  for  {y}  as 

{*•}  =  {TV<  (42) 

where  {y}  represents  the  amplitude  of  {y}. 

Substituting  Eq.  (42)  into  (41),  we  can  derive  the  characteristic  equation  for  the 


eigenvalue  A 

4 

where  Slj  are  given  by 

=  0 
i= o 

(43) 

(2<)  =  011022  “  012021 

(44a) 

(2l  =  <*11022  +  <*22011  —  <*12021  ~  <*21012 

(446) 

=  <*11<*22  “  <*12<*21  “  011  “  022 

(44c) 

SI3  =  —On  —  <*22 

(44  d) 

(24  =1- 

(44e) 
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The  stability  theorem  for  linear  systems  states  that  a  linear  system  is  stable  if 

and  only  if  the  roots  of  the  characteristic  equation  all  lie  in  the  left  half-plane,  excluding 

the  imaginary  axis.  To  determine  whether  the  polynomial  given  by  Eq.  (43)  has  all  its 

roots  in  the  left  half-plane  without  actually  solving  for  all  the  roots,  we  use  the  Routh 
22  23 

table,  ’  and  derive  four  parameters  for  the  polynomial  given  by  Eq.  (43).  For  the 
roots  of  this  polynomial  to  be  confined  in  the  left  half-plane,  excluding  the  imaginary  axis, 
these  four  parameters  must  all  be  strictly  positive.  Such  a  requirement  yields  the  following 
four  conditions 


Qo  >  0 

(45a) 

(2s  >  0 

(456) 

(22 (2s  —  (2i  >  0 

(45c) 

(2i((2j(2s  —  (2j)  —  (2q(2|  >  0. 

(45d) 
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Consequently,  the  linearized  system,  Eq.  (41),  is  stable  when  all  four  conditions  given 
by  Eq.  (45)  are  met.  As  an  example,  let  us  consider  the  first  two  conditions.  Substituting 
Eq.  (44)  into  (45)  yields 

$11*22  +  $22*11  -  $12*21  -  $21*12 


$11$22  -  $12$21 
X11X22  -  X12X21  >Q 


>0 


(46a) 

(466) 


$11$22  -  $12$21 

where  tyjj,  and  Xij  represent  the  mass,  damping,  and  stiffness  per  unit  area  of  the 
plate  and  are  defined  in  Eq.  (29).  In  particular,  $ij  includes  the  added  mass,  *,,  includes 
the  added  damping,  and  Xij  includes  the  added  stiffness  induced  by  acoustic  radiation.  If 
the  effects  of  fluid  loading  and  mean  flow  are  neglected,  then  $ij  =  *,,  =  0,  for  *  ^  j, 
so  the  inequality  in  Eq.  (46a)  is  automatically  satisfied.  Also,  since  the  stiffness  matrix 
is  predominately  a  diagonal  matrix,  the  inequality  in  Eq.  (46b)  always  holds  true.  Thus 
without  fluid  loading  and  mean  flow,  a  linear  system  is  always  stable.  Equation  (46)  thus 
shows  that  when  the  effects  of  fluid  loading  and  mean  flow  are  large  enough  to  change  one 
of  the  inequalities  given  by  Eq.  (46),  then  the  amplitude  of  the  plate  flexural  vibration 
may  grow  exponentially  without  a  bound. 

One  advantage  of  the  stability  analysis  discussed  above  is  that  it  allows  one  to  relate 
the  structural  instability  phenomenon  directly  to  the  physical  quantities  possessed  by  an 
elastic  structure.  In  a  following  paper,  we  will  use  this  stability  analysis  to  analyze  the 
mechanisms  of  structural  instabilities  under  various  fluid  loading  and  mean  flow  conditions. 
In  particular,  we  will  plot  stable  charts  so  that  we  know  when  and  why  an  elastic  structure 
may  become  unstable. 


IV.  EXCITATION  FORCING  FIELD 

In  the  preceding  section,  we  presented  a  general  stability  analysis  and,  in  particular, 
derived  stable  conditions  for  a  linearized  system.  In  this  section,  we  consider  the  steady- 
state  response  of  a  fin  plate  excited  by  a  turbulent  boundary  layer. 
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It  has  been  shown  that  the  plate  flexural  displacement  w(x,y,  t)  can  be  expanded 
in  terms  of  the  orthogonal  base  functions  {!Vr(x)}T[C(t)]{W’(y)}  along  the  x-  and  y-axes, 
respectively,  where  the  coefficients  [C(t)]  are  determined  by  Eq.  (28).  For  convenience  in 
the  derivation,  let  us  assume  that  (C(<)]  can  be  expressed  in  terms  of  a  solution  matrix  [£] 

[C(t))  =  [£][/(<)]  (47) 

where  [£]  is  independent  of  time  t  and  can  be  determined  numerically,  [/(<)]  represents 
the  forcing  excitation  -  due  to  a  turbulent  flow.  The  elements  of  each  of  the  sub-matrices 
[/'*(*)]  816  given  by  Eq.  (30). 

Since  the  turbulent  boundary  layer  is  random  in  nature,  one  must  rely  on  the  statis¬ 
tical  properties  of  the  turbulent  flow.  However,  in  many  applications  a  turbulent  boundary 
layer  can  be  assumed  to  be  stationary  in  time  and  homogeneous  in  space,  so  that  it  can  be 
described  by  a  space-time  cross-correlation  function  which  decays  with  spatial  and  tempo¬ 
ral  separations  and  convects  with  flow  at  velocity  U.  Based  on  the  mathematical  model 

4  5 

developed  by  Maestrello,  ’  we  can  express  the  ensemble  average  of  the  cross  correlation 
of  the  excitation  forcing  field  due  to  a  turbulent  boundary  layer  as 

(P7-(x,  K,0,  t)p’T(x\  y\ 0,  ('))  =  r/({,rj,0,r) 

g,  ^  f  \ 

where  superscript  *  indicates  a  complex  conjugation,  £  =  x  —  x',  r}  =  y  —  y\  r  =  t  —  t\ 
6  is  the  boundary  layer  thickness,  A\  =  0.044,  A2  =  0.075,  A3  =  —0.093,  A4  =  —0.025, 
Ki  =  0.0578,  Ki  =  0.243,  =  1.12,  KA  =  11.57,  a !  =  50/(UjRi),  here  VJIU  is  the 

equivalent  incompressible  Reynold’s  number,  03  =  0.26,  Ut  is  the  free  stream  velocity,  and 
(p2)  is  a  measure  of  the  mean  square  intensity  of  the  forcing  field. 
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rUrtiaii 'SHt 


The  corresponding  cross-power  spectrum  can  be  obtained  by  taking  a  Fourier  trans¬ 
formation  of  T /((,  77, 0,  r).  Evaluations  of  this  Fourier  transformation  can  be  facilitated  by 
the  residue  theory,  and  the  result  is 

P,(£,q, 0,w)  =  (?)(6/Ut)e-M/(a's)e-M(a*s)e-i("'u  £  Aje-iKi"s/u.  (49) 

;=i 

Next,  we  rewrite  the  ensemble  average  of  the  cross  correlation  of  the  forcing  field  T / 
with  pr  expanded  in  terms  of  the  base  functions 

T/«,  ■>, 0,r)  =  {IV(x)}r(/1{W'(x')},  (50) 

where  [F]  is  the  cross-correlation  matrix  defined  as 


w  =  m))mv))m»')v\r{t))T)-  (si) 

The  corresponding  cross-power  spectral  density  function  can  be  obtained  by  taking  a 
Former  transformation  of  Eq.  (50) 

=  {W(x)}T(i  /jF]e-VT)Wx')}.  (52) 

Equating  Eq.  (52)  to  (49)  and  using  the  orthogonality  properties  of  the  base  func¬ 
tions,  we  can  obtain  the  cross-power  spectral  density  matrix  whose  (m,n)th  element  is 


given  by 


[ujmYU-ury 


Wk(x)Wt{x')dx'dx.  (53) 


Once  this  is  done,  we  can  proceed  to  derive  a  cross-power  spectral  density  function 
for  the  plate  flexural  vibration  response.  First,  let  us  define  a  cross-correlation  function 
for  the  plate  flexural  displacement 

r  w((,rj,r)  =  (w(xty,t)w*(x',y',t')).  (54) 
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Second,  we  expand  to(x,  y,  t)  into  the  orthogonal  base  functions  using  Eq.  (26)  and  replace 
the  coefficient  matrix  [£?(<)]  by  the  solution  matrix  [£].  Doing  so  leads  to 


r-U.x.T)  =  W*)n£]W)TW*')} 


The  corresponding  cross-power  spectrum  of  the  plate  flexural  displacement  is  thus  found 
to  be 


{W(xo))TV,,u(U,r,,0^){W(zo))dx0dxa\ [£*f{WV)},  (56) 


where  P/,*i(£,»?,0,u>)  is  given  by  Eq.  (53). 

Equation  (56)  demonstrates  that  the  power  spectral  density  function  of  the  plate 

response  is  directly  related  to  the  power  spectral  density  function  of  the  excitation  forcing 

\ 

field  due  to  a  turbulent  boundary  layer.  The  auto-power  spectral  density  of  the  plate 
flexural  displacement  can  be  obtained  simply  by  setting  £  =  rj  =  0  in  Eq.  (56). 

V.  RADIATED  ACOUSTIC  PRESSURE 

In  this  section,  we  derive  a*  formulation  for  estimating  the  power  spectral  density 
function  of  the  radiated  acoustic  pressure  based  on  plate  flexural  vibration  responses.  First, 
let  us  define  a  cross-correlation  function  for  the  radiated  acoustic  pressures  measured  at 
two  different  points  in  the  fluid  medium 


r  =  <p(*,V*M)p*(*',v',  *',<')>,  (57) 

where  p*  is  the  complex  conjugate  of  p  given  by  Eq.  (14)  or  (17),  deper  ding  upon  whether 
the  flow  is  supersonic  or  subsonic.  Next,  we  replace  p  in  Eq.  (57)  by  its  Fourier  transfor¬ 
mation  pair 

(f *'»*')>=  /  /  ( p(*,V,z,w)pm(x\y'yz\u,))ci*e-™'t'dw<Li (58) 

t/— OO  J — oo 
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^•Jb':''5t, 


^A~’‘- 

_.,  ''-f‘^;r  . . \ 


■i ,+.■•■'>•*  -> *i '»«•'  /‘t*-  '.iii^L '  W~.".  '  rv  ~  nfi  r  '>n“"yn""Vi  rir tV'Sii r 


?.'>'(*  •••■  * 


The  frequency  domain  solution  for  the  radiated  acoustic  pressure  p  can  be  obtained  by 
substituting  Eq.  (7)  into  (13) 


#*,»,*,«)  =  ^jj(-iw  +  t r^jjo  l  G(x,y,z\x,,Vl,0) 

x  ^-iu  +  U u>(xi » Vi ,  w)dxi  dyi .  (59a) 


In  a  similar  manner,  we  can  write 


pV, »'.*>')  =  -(^j("  + 17^?)  jf  £ GV. y\ z’\xz,y2, 0) 

x  4-  17 tD*(x2,y2>w)dx2<iy2?  (59b) 

where  G*  is  the  complex  conjugate  of  G  defined  in  Eq.(9)  for  a  supersonic  flow  or  (11)  for 
a  subsonic  flow. 

Substituting  Eq.  (59)  into  (58)  and  rearranging  terms,  we  obtain 

2  r*  [L  r°° 

(p(x,y,x,t)p*(x',y',x',0)  =  j0  ill  J-ooJ-oo 

+ u  (*>'  +  u ^7)  G(x>  v>  *lxi  ’  y 1  >  °)G* (x<>  y' » z>  I**  >  V2  >  °) 

xe-*"te,w  *  dwdu)' dx\dxidy\dy%,  (60) 

where  the  triangle-bracketed  term  on  the  right  side  of  Eq.  (60)  can  be  rewritten  as 
(to(x  1 ,  yi ,  w)tD*  (xa ,  y2 ,  w' )) 

=  7^75-/  /  Mxi,yi,ti)u>*(x2,y2,f2))e",u,<1e'w  hdtidt9 

V2*)  «/ —00 «/ —00 

=  7 Tv  T  T 
(2*r  J-00  7-00 

=  (Sj /"  r.(£,,i)',r')e-^''<(u. -a,')dr',  (61) 
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where  £'  =  x\  —  x2,  r}'  =  y\  —  y2,  and  t'  =  t\  —  t2.  In  carrying  out  the  integration  with 
respect  to  t2  in  Eq.  (61),  we  have  used  the  definition  of  the  cross- correlation  function  rw 
(see  Eq.  (55))  and  the  property  of  the  Dirac  delta  function.^ 

Substituting  Eq.  (61)  into  (60)  and  using  the  definition  of  the  cross-power  spectrum 
for  Vw,  we  obtain 

Wx,„,m)pv,vv,<')> =<£777771 

+  +  ^^7)  G(x,y,x|xi,y,,0)G*(x',y',x'|x2,y2,0) 

where  ,r)'  ,w)  is  given  by  Eq.  (56). 

The  cross-power  spectral  density  function  for  the  radiated  acoustic  pressure  can  now 
be  obtained  by  taking  a  Fourier  transformation  of  Eq.  (62).  Using  the  property  of  the 
Dirac  delta  function,  we  can  carry  out  integrations  with  respect  to  t  and  u>  exactly  and 
obtain 

*-(£777 

+  U ^iu>'  +  U -9-^j  G(x,  y, z|x! , yj , 0)G*(x', y\ z' |x2, y2, 0) 
x  +  U (m  +  V Vw{? ,ri' ,u)dxidx2dyidy2.  (63) 

Equation  (63)  shows  that  the  cross-power  spectral  density  function  of  the  radiated 
acoustic  pressure  is  directly  related  to  the  cross-power  spectral  density  function  of  the 
plate  flexural  displacement.  If  we  substitute  Eq.  (56)  into  (63),  we  can  obtain  an  explicit 
formulation  for  V  in  terms  of  the  cross-power  spectral  density  function  of  turbulent  flow. 

It  is  emphasized  here  that  turbulent  flow  is  random  in  nature,  and  therefore  deter¬ 
ministic  formulations  for  calculating  the  resulting  structural  and  acoustic  responses  cannot 


Vw(t',ri',w)etu’rdwdxidx2dyidy2,  (62) 
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be  found.  The  formulation  given  by  Eq.  (63)  represents  a  new  way  of  estimating  the  radi¬ 
ated  acoustic  pressure  from  an  elastic  plate  excited  by  turbulent  flow.  In  deriving  Eq.  (63), 
the  effects  of  coupling  among  structural  modes,  radiated  acoustic  pressure,  and  turbulent 
flow  are  all  taken  into  account. 

VI.  CONCLUDING  REMARKS 

The  cross-power  spectral  density  function  of  the  radiated  acoustic  pressure  from  a 
finite  plate  subject  to  turbulent  flow  is  obtained  and  expressed  in  terms  of  the  cross-power 
spectral  density  function  of  the  turbulent  boundary  layer  in  an  explicit  form.  In  deriv¬ 
ing  this  formulation,  the  effects  of  stretching  due  to  in-plane  force  and  those  of  coupling 
between  structural  vibration  and  acoustic  radiation  are  taken  into  account.  A  general  sta¬ 
bility  analysis  is  given  by  using  the  basic  existence-uniqueness  theorem.  In  particular,  the 
stable  conditions  for  a  linearized  system  are  derived  by  using  the  Routh  algorithm.  Two 
sources  of  instabilities  are  found  to  be  attributable  to  the  added  damping  and  stiffness  due 
to  fluid  loading  and  flow.  The  effects  of  the  added  damping  and  stiffness  increase  with  the 
Mach  number  of  the  mean  flow,  the  former  increasing  linearly  and  the  latter  increasing 
quadratically. 
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